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Abstract 

The phase diagram of the quark-gluon matter evolution is presented for 
the SU(3)-model with a new phase of heavy deconfined quarks which ex- 
ists in a rather wide range of temperatures and densities. Fitting the chiral 
phase transition data to fix the model parameters we establish another (de- 
confinement) phase transition which separates a new phase from hadronic 
matter. The parameters and properties of the phase diagram are discussed 
in comparison with lattice and other results. 



^Research is partially supported by " Volkswagen-Stiftung" 

^Permanent address: Department of Theoretical Physics, P.N.Lebedev Physical 
Institute, Russian Academy of Sciences, 117924 Moscow, Russia. E-mail address: 
kalash@ td .lpi.ac.ru 



1 



At present for many applications which are based on hot QCD predic- 
tions and for future experiments it is very essential to establish a reliable 
(//, T)-phase diagram of the QCD-thermodynamical evolution and to cal- 
culate quantitatively the parameters which determine its phase properties. 
However this is not the case today since the properties of many (at least, 
two) phases are not elaborated in detail. Only the quark-gluon plasma phase 
(QGP-phase) which takes place at high temperatures and densities is reliably 
studied since due to a small coupling constant it can be treated perturba- 
tively. This is the phase of the "light" and deconfined quarks and ghions 
which are massless in the first approximation. We know as well that there is 
a phase which presents hadronic matter (H-phasc) where quarks and gluons 
are confined to form the observable nuclear matter. However, there is a ques- 
tion about the intermediate phase (Q-phase) which should present heavy and 
deconfined quarks. This phase is very natural and has a good physical mean- 
ing but till now it is not established firmly. Of course, it has been studied in 
many papers (in [1] and then in other papers) but no reliable predictions have 
been given. Today, the interest in this problem was revived again by paper 
[2] where the intermediate Q-phase is explicitly found and its parameters are 
determined within the bag model ideology. 

The goal of this paper is to derive the QCD-pattcrn of the Q-phase di- 
rectly from the Lagrangian approach using the standard calculations within 
the temperature Green function technique. Fitting the chiral phase transi- 
tion data to fix the model parameters we establish another (deconfinement) 
phase transition which separates the Q-phase from hadronic matter. Here 
we consider that the quark mass sharply arises after crossing the chiral phase 
transition line and the Q-phase demonstrates itself as a phase of heavy quarks 
with a rather strong interaction. The deconfinement phase transition occurs 
when the temperature decreases and within our scenario it does not coincide 
with the chiral transition. There are two well-separated phase transitions 
and below we specify their parameters. 

The QCD Lagrangian is usual and in covariant gauges it has the form 

>C= - \Glj' + Nfi,Yi,{d,-Ug\'^V^)+m]i, 

- i,Nfi^^,i^-r^{d,V^f + C\d,6'^' + gr"^V^)d,C' (1) 
where G^^ = d^V;^ - dyV^ + ^/"'"^V;^^; is the Yang-Mills field strength; 
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is a nonabelian gauge field; ■?/'(and ip) are the quark fields in the SU(N)- 
fundamental representation (|A" are its generators and /"^'^ are the SU(N)- 
structure constants) and C" (and (7") are the ghost Fermi fields. In Eq.(l) // 
and m are the quark chemical potential and bare quark mass, respectively, 
Nf is the number of quarks flavours and a is the gauge fixing parameter 
{a = 1 for the Feynman gauge). The metric is chosen to be Euclidean and 

The partition function is determined by the functional integral in the 
standard manner [3] 



Z^N j Viip, C, y)exp J d^x J dx4jC{ip, C, V) 



(2) 



where (3 is the inverse temperature and we calculate the thermodynamical 
potential for hot QCD following the formula 

n = -TlnZ. (3) 

The two-loop approximation for Q is considered at first but then we go beyond 
the perturbative theory. In the two-loop approximation the four diagrams 
(or three ones if the axial gauge is used) should be calculated. The calcula- 
tions are performed with the aid of the bare Green and vertex functions in 
the Matsubara technique. The Feynman gauge is more preferable although 
the final result docs not depend on the gauge choice. The renormalization 
prescription is well-known [4,5] and due to the quark loop (where is the mass 
renormalization) it is not trivial and more complicated then for hot gluo- 
dynamics. Such calculations were made many years ago in [4] for hot QED 
where the thermodynamical potential with nonperturbative corrections was 
found and then in [5,6,7] for hot QCD. In a more general case, when the 
quark mass and chemical potential m, 7^ 0, the two-loop thermodynamical 
potential was calculated in [7] and this expression is given by 

Q(^) _ n\N' - 1) g\N' - 1)N _ NNf_ J dp , 



V 45 144 37r2 J 

^ 

7V^(7V2 - 1)^2^2 Jd^ Nf{N' - 1)^^ Jdpda 

P ^ 
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X 



+ 



pq 



^ EpEq^ — m — pq 
EpE^ - m2 + pq^ 



/ - + , + , "^1 £^p^q + m' + pq 



pq -Bp-Bq + m? - pq^ 



(4) 



where rir 



n' 



rip and np 



[exp/9(i?p ± yU.) + 1] ^ are the usual quark 

2 I ^2\l/2 



occupation numbers. Here Ep = (p^ + m'^y^'^ is the quark energy. 

Our goal is to go beyond the perturbative expansion. There are several 
ways of doing it but phenomenologically (as a convenient fit) it is possible 
to replace g'^ by the running constant g'^{Q^) which further can be used for 
any temperatures and densities. Here we put our fit into agreement with the 
chiral phase transition data and then the Q-phase with the massive quarks 
will be investigated to find the deconfinement phase transition. The proposed 
fit has a rather simple form 



hi In 



20 



hi 



(5) 



1 + Q 16&0 



where feg = (HA^ - 2A^/)/3 and hi = {3AN^ - ISNNj + 3iV//iV)/3 are the 
standard renormalization group coefficients (for other more complicated fit 
see [8]). The ansatz (5) correctly reproduces the high momentum (temper- 
ature or density) behaviour and demonstrates the interaction strengthening 
when a small momentum region is considered. In this region where the in- 
frared divergencies are dominant the canonical dimensions arc removed from 
the theory (like in QED see e.g. in [3]) and anomalous ones arise. 

At first, we study the simplest case within the SU(3)-model when jj,,m — 
to check the fit (5) and to fix the parameter A. When the temperature de- 
creases the chiral phase transition should occur and namely this case is more 
convenient for lattice simulations and for analytical calculations as well: here 
all integrals in Eq.(4) are calculated exactly. The two-loop thermodynamical 
potential is given by 



V 



45 



■[{N^ - 1) + ^NNf] + g^T^ 



however Eq.(6) being a pure perturbative expansion does not demonstrate 
any phase transition [9]. The useful model arises only after replacing g"^ with 
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the aid of Eq.(5) and we study its properties below. The phase transition 
is determined by the condition {—p = Q{Qc) = 0) which here results in the 
equation 



6g In 1 + Qc 



■ bi , 
166^' 



144 [(iV2 



45{m - 1){N + iNf) 



(7) 



where Qc — TJA is a scaled variable and N — 3. Solving this equation one 
finds the critical Qc-points for each Nf as follows (see Table 1). 



SU{3) 


Nf = 2 


Nf = 3 


Nf = A 


niMeV) 


0.704 
186 


0.529 
160 


0.346 
130 



Table 1: The data summary for the chiral phase transition within the SU(3)- 
model. Here A = 222MeV to reproduce = 260MeV since Ql=1.375 for 
Nf — 0. For comparison, see lattice data in [10] and about A in [11]. 



Of course, the data in Table 1 can be improved but the order of their values 
and the tendency (the monotonic decreasing with Nf) are correct. The 
proposed model reliably reproduces the chiral transition temperature for any 
Nf. 

Now the model is fixed and we study the Q-phase with massive quarks to 
establish the dcconfinement phase transition when the temperature becomes 
lower. These calculations are more complicated since all integrals in Eq.(4) 
are not treated exactly when the quark mass m is non-zero. However for 
fjL — O Eq.(4) keeps a rather convenient form 

V 

+ 
+ 

since all integrals in this case depend on the scaled parameter uj = m/T. 



ttHN^-I) ^, _ 2NNfT^^^ 



45 



37r2 



g\N'-l)N^,^Nf{N'-l)g'T^ 



UA 



12 



UJ 



^ ( 3 liiu) + hsiu] 



(8) 
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These integrals are 



/2M 







Vx'^ + 1 exp(a;Va^^ + 1) + 1 
x'^dx 1 



\/a?~+Texp(c<Jv^a?~+T) + 1 





7 xdx ydy (,x - 1/)^ 

= / ^=^^^=^^in- 







X I I (9) 

exp(a;\/a?"+I) + 1 exp(a;\/|/^ + 1) + 1 

and they have to be calculated numerically. But, there is a problem with 
the Js-integral which is not prepared properly for numerical calculations. 
However this integral being rather small for uj >\ does not affect essentially 
the final result and due to this fact its logarithm can be shghtly redefined 

^ + (10) 

{x + yY (x + y)2 + 0.01 ^ ^ 

to make this integral convenient for the numerical calculations as well. 

Let us find the critical point which determines the deconfinement phase 
transition in our scenario. We consider the SU(3)-model with Nf={2 and 
3): the experimental data seem to be between these values. It is also as- 
sumed that all constituent quarks have the same masses which are gener- 
ated due to the chiral phase transition and have the magnitude of ISOMeV 
< m <350MeV. This quantity being a free parameter is fitted with the aid 
of the a;-variable to define another (deconfinement) phase transition. The 
equation which determines the critical temperature of this transition is es- 
tablished through the "standard" requirement (— p = ^{Qd) — 0) and using 
Eq.(8) its explicit form is given by 



45 ' TT^ 7r2 




( ^ + ^ [ 3 I2M + jhi^d) ] ) \ (11) 
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where Ud = ni/Td and Qc = T^/A. In Eq.(ll) all integrals are calculated 
numerically for the chosen cu^- values and then the equation obtained is solved 
to find the Q^-parameter. The results are summarized in Table 2. 



SU(3) 


Nf = 2 


Nf = 3 


Qi 


Td{MeV) 


m{MeV) 


Qi 


Td{MeV) 


m{MeV) 


0.7 


0.577 


169 


118 


0.400 


140 


98 


1.0 


0.540 


163 


163 


0.361 


133 


133 


1.3 


0.523 


160 


208 


0.339 


129 


168 


1.9 


0.543 


164 


312 


0.343 


130 


247 


2.5 


0.616 


174 


435 


0.395 


140 


350 



Table 2: The data summary for the dcconfincment phase transition within 
the SU(3)-model for /i = 0. We keep A = 222MeV to estimate Ta. For 
comparison see [2] where the similar results are discussed. 



Prom Table 2 we see that indeed Td=166MeV (for Nf^2) and T<i=140MeV 
(for Nf—3) seems to be the upper limit for this value since the quark mass 
is of the order of 350Mev or less. But in any case Td 7^ Tc (for Tc see Table 
1) and this fact is more important. These results completely confirm the 
scenario where the Q-phase separates QGP from the usual hadronic mater: 
the two well-separated phase transitions are established with AT^nax — Tc — 
Td ~30MeV. 

Below the case T = is investigated. Treating this case we start again 
from the two-loop approximation for Q which for T = is simpler since all 
integrals can be calculated analytically. Here one finds that the two-loop 
thermodynamical potential for cold quark-gluon matter has the form [7] 



V 



+ 9 



127r2 
MN' - IK 



2 2 u 



i + y/r^ 

6 



647r4 

- 2{1-9Y} (12) 
where 9 — m/ /i is a new parameter and the link between /i and the quark 
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density p is given by 



(13) 



At first we consider the chiral phase transition which occurs in the quark- 
gluon plasma with the hght (here massless) quarks. In this hmit Eq.(12) 
becomes very simple and, for A?" = 3, it is reduced to [6] 



V 



47r2 



1 - 



(14) 



where Q = h/hq. The phase transition occurs when p = and the equation 
which results from this condition has the form 



2 b?, In 



l + Qc 



2(4-) 



(15) 



where all notations are the same as above. This equation is solved numerically 
to find Qc and the results are summarized in Table 3. 



SU{3) 


Nf = 2 


Nf = 3 




0.232 
700 


0.195 
618 



Table 3: The data summary for the chiral phase transition within the SU(3)- 
model when T = 0. Here we choose /xq = 1.4GeV to estimate the quark 
chemical potential Pc in accordance with [2,12]. 

The deconfinement phase transition occurs when the quark density decreases 
and quarks become massive: their mass is generated on the chiral phase tran- 
sition line. The critical density is determined in the same manner through 
equation (-p = ^(Ql) = 0) whose explicit form is found by using Eq.(12). It 
is given by 



bi 



2 yr^(l - l^i) + 3 9j ]n- 



1 + 



bl In 



1 + 



2(-^) 



1 + 



1-9', 



(16) 



- 2(1 - ej) 



2\2 
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where Od = m/fid- When solving Eq.(16) one finds, at once, that there is an 
upper bound for 6'd-parameter: only 9d < 0.4 is acceptable to keep the right 
side of Eq.(16) positive. This bound means that the quark mass generation 
is suppressed when T = and we find ~ 60 — 80Mcv only. Of course, 
this estimate depends on the scale fiQ which we fixed above but in any case 
the chiral phase transition demonstrates itself in the (T = 0,/i)-rcgion very 
smoothly. Two values for the 6'd-parameter are considered below and for both 
cases we find approximately the same quark mass but completely different 
values of A// = A*c — (see the results in Table 4) . 



SU(3) 


Nf = 2 


Nf = 3 




Q'i 




mq{MeV) 


Qi 


^id{MeV) 


mq{MeV) 


1/8 


0.148 


538 


67 


0.106 


456 


57 


1/6 


0.098 


438 


73 


0.068 


365 


61 



Table 4: The data summary for the deconfinement phase transition within 
the SU(3)-model when T = 0. Here we keep /xq = 1.4GeV to estimate /id- 
For comparison see [2] where the similar results are discussed. 



From Table 4 one can see that the proposed scenario again presents two well- 
separated phase transitions with a rather broad Q-phase but, unfortunately, 
the mass generation in the case T = is not expressed sharply. 

To conclude two well-separated phase transitions have been established 
within our scenario. We find that AT ~20MeV on the axis /i = and 
the larger range A//, ~( 160-260) MeV on the axis T = 0. The intermediate 
(/X, T)-region is not investigated but it is obvious that the scenario should 
be close. Of course, the present scenario is only a model which should be 
improved in a selfconsistent manner but it being in agreement with the chiral 
phase transition data is a solid basis for the theory to come. For removing 
any arbitrariness it is necessary to solve nonperturbatively the Schwinger- 
Dyson equation for niq (e.g. see [13]) and to combine this result with the 
present scenario. On this way one determines both the dynamical quark mass 
mq(T,ix) and the deconfinement phase transition more exactly. The fit (5) 
which we used for the running coupling constant is also should be improved. 
But this question (the same as one about the higher order corrections) being 
very important and complicated requires additional investigations. Here we 
choose the simplest fit for g'^{Q'^) to perform the lowest order calculations. 
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